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We introduce a family of compatible Poisson brackets on the space of 2 x 2 polynomial matrices, 
which contains the Sklyanin bracket, and use it to derive a multi-Hamiltonian structure for a set 
of integrable systems that includes XXX Heisenberg magnet, the open and periodic Toda lattices, 
the discrete self-trapping model and the Goryachev-Chaplygin gyrostat. 

1 Introduction. 

The ingenious discovery of Magri [6j [7j that integrable Hamiltonian systems usually prove 
to be bi-Hamiltonian, and vice versa, leads us to the following fundamental problem: given 
a dynamical system which is Hamiltonian with respect to a Poisson bracket {.,.}o, how 
to find another Poisson bracket { , -}i compatible with initial bracket and such that our 
system is Hamiltonian with respect to both brackets. This, along with the related problem 
of classification of compatible Poisson structures, is nowadays a subject of intense research, 
see e.g. |6j U\ El EH] and references therein. 

In this paper we study a class of finite-dimensional Liouville integrable systems described 
by the representations of the quadratic r-matrix Poisson algebra, or the Sklyanin algebra: 

{T{X),T(fj,)} = [r(A - n),T(X)T(n) ] , (1.1) 

1 2 

HereT(A) = T(A) ® I , T(fj) = I® T(n) and r(A — //) is a classical r-matrix [8]-[TT]. 

The main result of the present paper is a family of the Poisson brackets {., which 
is compatible with the Sklyanin bracket (11.11) . in the simplest case of the 4x4 rational 
r-matrix 



r(\- fi ) = -*—IL, n 
A — fi 



( 1 \ 

10 

10 

\0 1/ 



rjeC, (1.2) 



and 2x2 matrix T(X), which depends polynomially on the parameter A 

™ " ( cIa] ) CM) 

a\ n + A 1 X n ~ 1 + ... + A n f]X n + B^- 1 + . . . + B n 
7 A n + dA™- 1 + ... + C n SX n + DiX n ~ l + ... + D n 
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The leading coefficients a, (3, 7, 5 and 2n coefficients of the det T(A) 

d(X) = det T(A) = (a«J - f3 1 )X 2n + QiA 2 ™" 1 + ■ ■ • + Q 2n . (1.4) 

are Casimirs of the bracket (II. II) . Therefore, we have a 4n-dimensional space of the coeffi- 
cients A{, Bi, Ci and Di with 2n Casimir operators Qi, leaving us with n degrees of freedom. 
For so-called open lattices independent Poisson involutive integrals of motion H° = Aj, 

1 = 1, . . . ,n, are given by the coefficients of the entry A(X): 

A(X) = a\ n + iJ°A n_1 H H°, {H°,H°} = 0. (1.5) 

In generic case integrals of motion are given by the coefficients of the trT(A): 

trT(A) = (a + 5)X n + H^ 1 + ■ ■ ■ H n , {H i ,H j } = 0. (1.6) 

These integrals of motion define two Liouville integrable systems, which are our generic 
models for the whole paper. Bi-hamiltonian description of these models gives rise to the 
bi-hamiltonian description of the Goryachev-Chaplygin gyrostat [5], open and periodic Toda 
lattice [9], inhomogeneous Heisenberg magnet [11] and the discrete self-trapping (DST) model 
0. 

2 The compatible bracket 

In this section, we describe the Poisson bracket compatible with the Sklyanin bracket. The 
Poisson brackets {., .}o and {., .}i are compatible if every linear combination of them is still 
a Poisson bracket. The corresponding compatible Poisson tensors Pq and Pi satisfy to the 
following equations 

[Po,P O ] = [Po,Pi] = [Pi,Pi] = 0, (2.1) 

where [., .] is the Schouten bracket [21 [6j [7]. Remind that on a smooth finite-dimensional 
manifold the Schouten bracket of two bivectors X and Y is an antisymmetric contravari- 
ant tensor of rank three and its components in local coordinates z m read 

[X, Yf k = -EUV+r fc ^ + cycle(z, j, k) . 
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2.1 Open lattices 

The Sklyanin bracket (jl.ip amounts to having the following Poisson brackets between the 
entries A(X), B(X), C(A) and D(X) of the matrix T(A): 

{A(X),A(fi)} = {B(X),B(fi)} = {(7(A), C(/z)} = {D(A), £>(//) } = 0, 
{B(X),A(n)} = ^L-[B{X)A{^) - B(riA{\j), 

{C(X),A(fi)} = Y ^-{c{X)A^) - C(v)A(X)), 

{B(X),C(n)} = (p(\)AQi) - D(vl)A(\)). (22) 
{B(A),D(/i)} = ^=5_(fl(A)D0*) - B(riD(\j), 
{C(X),D(»)} = j?L-(c{\)DQi) - Cfa)D(\j), 

{A(X),D(n)} = ^_(c(A)fl(M) - COi)B(A)), 

In fll.ip matrix r(A — //) satisfies the Yang-Baxter equation, which ensures the Jacobi identity 
for the brackets (12.21) . 

Proposition 1 The Sklyanin bracket U.1\) . I[K^) is compatible with the following bracket 
{•••I:-- 

{A(A), Ain)}, = {B(X), B(n)} x = {C(A), C(a*)}i = 0, 

{5(A), AOOh = -^(A5(A)A(/i) - ^)A(X)) - ^(A)A(/i), 
A — ft V J a 

{C(X), AMh = T ^-(xC(X)A^)-fiC^)A(X)) + H2 A(< x)A(fi), 
A — /i V J a 

{5(A), = JL-(\D(\)A(n) - nD(fx)A(X)) - IIa(X)A^), 

X — (J,\ J a (2.3) 

{5(A), £■(/,)}, = ^(B(X)D(f,)-B^)D(X))+ V A(X) (?-D{p) - -B(ji)) , 
X — J \at a J 

{C(A), £)(/*)>! = ^-(C(\)DM -C(n)D{\j) -r,A(X) (%(/.) - -C(jjl)) , 

{A(A), £)(/*)>! = Y^j^(c(X)B(fj,) - C(fi)B(X) S j - r)A(X) flflfo) - ^C(^) 
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and 



{D(\), D(fj,)} 1 = ^ ^D(A)-B(/f) — D(fj,)B(X)j — ^D(A)C(/^) - D{^)C{X)^ 

+ ^(B(X)C{fi)-B(ji)C(X)). (2.4) 



Proof: It is sufficient to check the statement on an open dense subset of the Sklyanin 
algebra defined by the assumption that A(X) and B(X) are co-prime and all roots of A(X) 
are distinct. 

This assumption allows us to construct a separation representation for the Sklyanin 
algebra fll.ll) . In this special representation one has n pairs of Darboux variables, A;, 
i — 1, . . . , n, having the standard Poisson brackets, 

{Xi, Xj} = {fit, Hj} = 0, {Xi, fij} = Sy, (2.5) 

with the A-variables being n zeros of the polynomial A(X) and the /i-variables being values 
of the polynomial B(X) at those zeros, 

A{Xi) = 0, & = r]- 1 In B{Xi), i = l,...,n. (2.6) 

The interpolation data (I2.6P plus n identities 

J B(A,)C(A i ) = -d(Xi) 

allow us to construct the needed separation representation for the whole algebra: 

A (X) = a(X-X 1 )(X-X 2 )---(X-X n ), 

B{X) = A(X) £ + ELi 



a (A — Aj)A'(Aj) 

d(A) + £(A)C(A) 

The coefficients of the determinant d(A) fll.4j) are Casimir elements for the both brackets 
{., .}o and {., .}i and, therefore, we can easy calculate the bracket {., .}i (I2.3p - (I2.4I) in (A, /i)- 
variables 

{•\, Xj} 1 = {/ij, = 0, {Aj, = Aj(5jj, (2.8) 

In order to complete the proof we have to check that brackets (12.81) is compatible with the 
canonical brackets (12. 5p . The compatibility of the brackets (I2.5p . (l2.8p implies the compati- 
bility of the brackets (I2.2p . fl2.3p and vice versa. 

The (A, /i)-variables (12.61) are so-called special Darboux-Nijenhuis coordinates [61 [7J [2] 
because 

. l\ p-( diag(Ai,...,A„) 

" TOy' 1 1 -diag(A 1 ,...,A n ) 



and the corresponding recursion operator N takes the diagonal form 



n ( d d \ 

iV = P 1 P -i = £A,(^dA, + -<a*,) 

i=l x ' 



(2.9) 



These Poisson tensors Po an d Pi satisfy to the equations (12.11) and the Nijenhuis torsion of 
N vanishes as a consequence of the compatibility between Po and Pi- 

Proposition 2 Brackets \2.5\) and Ii2.8\) between (X, p)-variables belong to a whole family 
of compatible Poisson brackets {., .} k associated with the Poisson tensors 



Ffe " iVfcFo - ( -diagCA},...,**) 



diag(Aj, . . . , A*) 




k = 0, 



In the matrix form, these brackets are equal to 



n. 



r(A),r(M) 



rf 2 ] (A, (j)T(X)T([m) -T(A)T( M ) rg ] (A, p) 



(2.10) 



Here 



/ 



1 o 



V 
\—fl 



1- 



M fc o 
A fc l- x 



,0 1 



and 



/ pg 1 oo 



o 







[fc] [fc] A fc - M fe „ 

x Pc Pb ~^r- 
\ 0, 



[fc] _ A fc X(A) p k X(p) 
Px 



s l S(X,p) = Us^(X,p)U. 



(2.11) 



w/iere X = B,C, D, 



A(X) A{p) 

is a difference of two polynomials, which are quotients of polynomials in variables X and p 
over a field. 



Proof: At k = one has p B = 0, p c = and p D = 0, so the bracket (12.101) coincides with 
the Sklyanin bracket (11. ip . At k = 1 we have 



[l] _ g(A - p) 



a 



fi] _ 7(A ~ M) 



a 



m _ gA - /i) 



and bracket (I2.10p coincides with the bracket (12. 3p . 

At k > 1 one can easily check that fc-th brackets ( 12. lOf) between polynomials A(X), B(X), 
C(X) and D(X) (EZJ) imply the brackets 



{Xi, Xj} k = {p l7 pj} k = 0, {Xi, pj} k = X k 5^ 



(2.12) 
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and vice versa. This completes the proof. 

To proceed further we need to recall that the normalized traces of the powers of N 

1 n 

J m = trace iV m = > A" 1 , m = l,...,n. (2.13) 

2m 

are integrals of motion satisfying Lenard-Magri recurrent relations [HI [7] 

iWi = 0, Xj i =P dJ i = P 1 dJ i - 1 , PidJ n = 0. (2.14) 
By definition ( 12 .TP polynomial 

n 

A(A) = aA n + A 1 X n - 1 + . . . + A n = a JJ(A - \) 

i=i 

is directly proportional to the minimal characteristic polynomial of iV (I2.9p 

n 

A^(A) = (det(JV - AI)) 1/2 = - \). 

i=l 

Since Hamiltonians H° (11.51) are related with integrals of motion J m (12 . 13[) by the triangular 
Newton formulas 

aJ 1 = H° l , a J 2 = H° + ^-, a J 3 = H° + H%Hf + . . . . 

As a consequence of the recursion relations (12.141) . the Hamiltonians H°, i — 1, . . . , n, satisfy 
the Frobenius recursion relations 

N* dH° = dH° l+1 - a- x A % dHl, , (2. 15) 

where iV* = Pq 1 P\ and H° +1 = 0. Such as Ai = H° a straightforward computation shows 
that they are equivalent 

N* dA(X) — A dA(X) + A{\)dA 1 . 

The special Darboux-Nijenhuis coordinates A,, fa are variables of separation of the action- 
angle type [2], i.e. the corresponding separated equations are trivial 

W,A J } = {J,,A i } = 0, i,j = l,...,n. 

We can introduce another separated coordinates Ui,Vi, which are the so-called Sklyanin 
variables defined by 

B(ui)=0, Vi = In A(ui) , i = l,...,n. 
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The separation representation of the algebra in (u, v)-variables has the form 

B(X) = /3(X-u 1 )(X-u 2 )---(X-u n ), 



A(X) = B(X)(^- + J2 



-VVi 



5 A d( Ui ) e" Vi 



C(A) 



A(X)D(X) - d(X) 



B(X) 

Substituting matrix T(X) (11.31) with these entries into the brackets {.,.}& (12.10!) at A; = 0, 1 
one gets that Ui, Vj coordinates are Darboux variables with respect to the Sklyanin bracket 

{ui, uj} = {vi, Vj} = 0, {ui, Vj} = Sij, (2.16) 

whereas the second brackets look like 

r e P-MVi) , , A'iuA A'tuA 

The corresponding separated equations 

n—l > 

{A(X),Uj} k = X k A(uj) 11^^' J = !,...,«. (2-1T) 

are linearized by the Abel transformation on the algebraic curve defined by e~ VVi = A(ui), 
see [91 131 H2] and references within. 

The special Darboux-Nijenhuis coordinates are dual to the Sklyanin variables. Namely, 
Aj, jUj are roots of polynomial A(X) and values of polynomial 5(A) at A = Aj, while Ui, Vi are 
roots of polynomial B(X) and values of polynomial A(X) at A = itj. 

2.2 Generic model 

There are many other Poisson brackets compatible with the standard one (12.51) . The main 
property of the proposed above bracket {., .}i (I2.3p - (I2.4I) is that 

{A(X),A( f i)} = {A(X),A( f i)} 1 = 0. 

It ensures that integrals of motion H° for the open lattices are in bi-involution 

{^;} = K>^}i = o - 

In this subsection we are looking for bracket {., .}'±, which has to guarantee the similar 
property for generic integrals of motion Hi (II. 6p from trT(A) 

{Hi,Hj} Q = {H^HjYi = 0, i, j = 1, . . .,n. 

Remind that {., -}o is the Sklyanin bracket (11.11) . which already has the necessary property 

{trT(A),trT(/x)} = 0. 

The following propositions can be ascertained by means of direct calculations. 
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Proposition 3 If a = 5 and f3 — 7 = in T(A) M.'d\) . then 

{trT(A),trT(/i)} 1 = 0. (2.18) 

So, the desired bracket {., .}[ may be obtained from the bracket {., .}i (I2.3I) - (I2.4I) by using 
special canonical transformations, which are generated by the suitable transformations of 
the matrix T(A). 

Proposition 4 The Sklyanin bracket U.l\) is invariant with respect to transformation 

T(A)^^T(A)F 2) y <=(^f)' <Xi,Pi,H,SieR, (2.19) 
where Vi^ are numerical matrices. If 

/3i72 + = 0, 

the bracket lji2.3\) - (2~4\ ) after transformation $2.19\) has the necessary property 

{trT(A),trT(/i)}i = 0. (2.20) 

We present an explicit form of the bracket {., .}[ in the Section 4 devoted to the periodic 
To da lattice. 



3 The Heisenberg magnet 

Another important representation of the quadratic algebra with the generators Ai,Bi, G{ 
and Di comes as a consequence of the co-multiplication property of the Sklyanin algebra 
(11. ip . Essentially, it means that the matrix T(A) (11. 3p can be factorized into a product of 
elementary matrices, each containing only one degree of freedom. In this picture, our main 
model turns out to be an n-site Heisenberg magnet, which is an integrable lattice of n sl(2) 
spins with nearest neighbor interaction. 

In the lattice representation the matrix T(A) (11.31) acquires the following form: 

T(A) = Lx(A - a) L 2 (\ - c 2 ) • • • L n {\ - c n ) , (3.1) 

with 

MA) = [ ( m) _ ? fl ( ra) \ + s (A y m =!,..., n. (3.2) 

Here are dynamical variables, c m are arbitrary numbers and i = \J — 1 . 

Substituting matrix (13. ip into the Sklyanin bracket (11. ip and brackets (I2.3p - (I2.4I) at 77 = i 
one gets canonical brackets on the direct sum of sl(2) 

{^\sf)=e m st\ (3.3) 
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and second compatible brackets 

K".*? }, = •?>(*. - 4 m) ). {r ,>? } l = (^ m,) ) y , m*t, 

where is the totally skew-symmetric tensor and 

/ i^H„W , „M «)\ • (m) tt) (m) U) i« (m) r«i w \ 

H 6 3 6 3 2 *2 J 1A 2 6 1 6 3 6 3 1A 3 ^1 lb 2 ) 



pipit) 



1 V*3 6 3 1 5 1 ,1 *3 ^ A l 1,b 2 J 



is 1 s 2 + s 3 s 3 



V -i(4 m) + i4 m) )^ (4 m) + i4 w) )4^ -i(4 m) + i4 m) ) - J 

The corresponding Poisson tensors Po and Pi are degenerate and, therefore, the Hamiltonians 
H° satisfy the Frobenius recurrence relations (12.151) in the following form 

PidP° = P (dH° +1 - AidHl) , i = l,...,n, (3.4) 

where H° +1 = and A$ = H° are coefficients of the polynomial A(X). The first integrals of 
motion are 

a? = X> - 4 m) ) ■ m = E(.f» - i4 m) ) (4" + i4") - 5 X>» - 4"*) 2 + ^ ■ 

m=l m>£ m=l 

Such as a = 5 and = 7 = we can use these brackets for the open and periodic lattices 
simultaneously. It means that Hamiltonians Hi (I1.6P from the trT(A) satisfy the Frobenius 
equations (13 .4p too. 

4 The Toda lattices 

The Toda lattices appear as a specialization of our basic model when the parameters are 
fixed as follows: 

= 7 = 5 = Q and det T{\) = 1. (4.1) 

We also put a = 1 and rj = — 1. In the lattice representation, the monodromy matrix T 
(I1.3P acquires the form 

T(A) = Li(A) • ■ ■ L n _!(A) L n (A) , U = ( X ~ q Pl ~f ) . (4.2) 

Here p^, are dynamical variables. 

4.1 Open lattice 

Substituting matrix T(A) (14.21) into the brackets {., .}k (12.101) at k — 0, 1 one gets that the 
Poisson tensors Po and Pi in (p, q) variables take the form 

^0 = £#- a /T' (4 - 3) 

Tl 1 r\ r\ Tl r\ r\ Tl r\ r\ 

P A A A 
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Namely this bi-hamiltonian structure of the open Toda lattice was obtained in pQ. 

For the open Toda lattice the Hamiltonians H° from the A(\) = X n + i^A n_1 + . . . + H° 
satisfy the Frobenius relations (I2.15p . The first integrals of motion are equal to 



n i n n—1 _. / n \ ■ 

8=1 i=l i=l \i=l / 



(4.4) 



The Sklyanin variables iti, fj are introduced as before: 

B( Ui ) = 0, v i = -ij- 1 \nA(ui), z = l,...,n-l, (4.5) 

the only difference now is that this gives only n — 1 instead of n separation pairs. The missing 
pair of canonical variables is defined as follows: 

n 

v n = In &i = -q n , Un = -a 1 = ^2p i . (4.6) 

i=i 

The separation representation of the algebra in (u, i>)-variables may be found in [9], [13]. It 
is easy to prove [13] that (it, u)-variables are Darboux variables 

« 

lo = Pq 1 = dui A dvi, 
i=i 

and the only nonzero second Poisson brackets are 

{uj,Vi} 1 = u i Sij, {u n ,Ui} 1 — —e " wjui) > { u m Vi} 1 = — e " bu^h) ' 

n 

{Vn,Vi}i = -1, {u n ,V n } 1 = ~22 U i- 

i=l 

Remark 1 From the factorization (14. 2 j) of the monodromy matrix T(A) one gets 
B n (A) = -e*»4i-i(A) = -e**4i-i(Ai) = 0. 

This implies that for the (n—1) -particle chain special Darboux-Nijenhuis variables Xj coincide 
with the Sklyanin variables Uj, i = 1, . . . , n — 1 for the n-particle chain. 

4.2 Periodic lattice 

For the Toda lattice a ^ 5 and, therefore, in order to get new bracket {., .}^ with the the 
necessary property (12. 20j) we have to apply transformation (I2.19P to the initial bracket {., - }i 
fl2T3> fl2~l. If we put 



I 1 )' and V *=[ \ \ 
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(4.7) 



then one gets the following brackets between the entries of T(X): 

{A(X), A(n)}[ = rj(B(X)C(fi) - Bfr)C{\)) , {5(A), 5(//)}; = 0, 

{A(X),D(n)}[ = -^-(C(A)BOi) - C(ji)B(X)) 

A — fi 

{B(X), Bin)}', = V (B{\)D(n) - B(jj)D(X)), 
{C(X),C^)}[ = rj(C(X)D(fi) - C(jj)D(X)), 

{D(X), B0*)K = J!t-(B(X)D(ji) - B(fi)D(X)) 

A — jl 

{D(X), C ([!)}[ = f^(C(A)5(/i) - C(jj)D(X)) 

A — /i 

and 

{A(X),B^ yi = --^^ + (5(A) - C(X))BM 

{A(X)M»)Yi = ,( ~ (A)CM - V(C(X)D(,) + (5(A) - C(X))C(,)), (4-8) 
{5(A), C( M )K = ^^(>)-M(A) g (,)) _ ^ (B(A)2?(/i) _ D(A)c(At) + D(A)D(/i)) . 

According to proposition 4 these brackets have the necessary property (12.20)) . 

Substituting matrix T(X) (14. 2 p into the brackets (14. 71) - (14. 8 p one gets that the Poisson 
tensor P[ in (p, q) variables takes the form 

Yl 1 r\ r\ TL r\ r\ TL r\ r\ 

p, qi-q i+ i O O O O ST^ O O 

(e qi ® A ^ + e 9n ® A ^ ^ 
^ V % dpi dqi dp n ) 

For the periodic Toda lattice the Hamiltonians Hi and H 2 from the trT(A) = X n + HiX n ~ l + 
. . . + 5 are equal to 

n i ra n i / n \ ^ 

Hi = H{ = -Y^Pi, H 2 = H° + = - 5> 2 + - - ( J> J , (4.9) 

i=l i=l i=l \i=l / 

where = These Hamiltonians Hi, i — 1, . . . ,n, form the Frobenius chain 

7V*cL5; = dH l+1 + c t dH x , with H n+l = . (4.10) 
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Here iV* = P P[ and q are coefficients of the minimal characteristic polynomial of the 
recursion operator 

Ajv(A)= (det(iV-AI)) 12 = A n -(c 1 A n - 1 + ... + c ft ) J (4.11) 

which can be defined directly via the entries of the matrix T(A) 

Ajv(A) = A(\) + B{\) - C(A) - D{\) . (4.12) 

Remark 2 Transformations (I2.19P of the matrix T(A) give rise to canonical transformations 
in the phase space. As sequence tensor P[ (14.91) coincides with tensor P\ (14. 3 p after the 
following canonical transformation 

Pi ^ Pi + e~ gi , p n — ► p n + e qn , 

which identifies coefficients q with integrals of motion for the open Toda lattice q = —H°. 

5 Integrable DST model 

The integrable case of the DST (discrete self-trapping) model with n degrees of freedom was 
studied in [5]. It appears as a specialization of our basic model when several parameters 
vanish: 

f3 = 7 = d = and Q s = 0, j = l,...,n-l. (5.1) 

We also put a = 1 and r] = —1. In the lattice representation, the matrix T(A) (11.31) acquires 
the form 

T(A) = Li(A - ci) L 2 (A - c 2 ) • • • L n (A - c n ), with L*(A) = ( A ^ . (5.2) 

Here Pi,qi are dynamical variables, whereas b and are numbers entering into the Casimir 
function (Oj) 

d(A) = detT(A) = 6 n (A - ci)(A - c 2 ) • • • (A - c n ). 

Substituting matrix T(A) (15.21) into the Sklyanin bracket fll.ip and into the brackets (14 . 7ft - 
(14. 8 p one gets canonical brackets 

{Pi>??}o = <y v» {ft>9j} = {Pi'PAo = °' i,J = l,...,7l. 

and quadratic brackets 

{?i,9i}i = -QtQj, {Vi,Pj}i = Qi p j - c i^i^ i > J 

{Pi , Pi) i = -PiPj , {Pi , = QiPj - b 5 i+lj , 
where Qx — q± + 1, P n — p n + b and Qi = qi, Pi = pi for other values of index i. 
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As above the Hamiltonians ifj, i — 1, . . . , n, from the trT(A) = A ra + Hi\ n 1 + . . . + H n 
satisfy the Frobenius relations (I4.10p . The two first Hamiltonians of the system are 

n 

Hi = -^(giPi- Cj), (5.3) 
i=i 

n 

H 2 = ^(QiPi ~ Ci)(qjPj - cj) ~b^2 liPi+i, Pn+i = Pi- 

i>j i=l 

The Sklyanin variables i — 1, . . . ,n, are introduced by the same formulae as for the 

Toda lattice, cf. flU} and 



6 The Goryachev-Chaplygin gyrostat 

Let us consider the matrix T(X) introduced in [I] 



T(A) 



A 2 - 2A J 3 - J\ - 
(xi - \x 2 )\ - X 3 ( Ji 



J 2 



(x 1 + ix 2 )A - x 3 ( Ji + iJ 2 ) 



iJo 



-x\ 



(6.1) 



Substituting matrix (16.11) into the Sklyanin bracket (11 .11) and brackets (I2.3p - (I2.4I) at r) — 2i 
one gets canonical Poisson tensor on the dual space of Euclidean algebra e(3) 



(6.2) 



(0 











x 3 


-x 2 \ 


* 








-X 3 





Xi 


* 


* 







—X\ 





* 


* 


* 





h 


-h 


* 


* 


* 


* 





Ji 


V* 


* 


* 


* 


* 


o / 



and the following quadratic tensor 



Pi 







X 3 X 2 


-X2J1 


-X2J2 


x 3 J 2 - 2x 2 J 3 \ 


* 





-X3X1 


XiJ\ 


X1J2 


2x 1 J 3 - X3J1 


* 


* 











-x 1 J 2 + X 2 Jl 


* 


* 


* 





_ 72 _ 72 


-J3J2 


* 


* 


* 


* 





J1J3 


V* 


* 


* 


* 


* 


/ 



(6.3) 



These tensors satisfy equations (12.11) at any values of the Casimir functions 

Ci = x\ + x\ + x 2 , C 2 = X1J1 + X2J2 + x 3 J 3 . 

However, in the proposed method coefficients of detT = —C\\ 2 + X3C2A have to be the 
Casimir functions and, therefore, we have to put C2 = 0. As sequence, we have {^4(A), A{ji)}i = 
at C2 = only. 
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Remark 3 Solving equations P^dH^ = (Pi + o.Pq)H° at arbitrary values of C Xj2 one gets 

H? = J 3 , H° = Jl + Jl + 2Jl + aJz. 

Here H 1 ^ is a kinetic part of the Hamiltonian for the Kowalevski gyrostat, which may be 
studied by using 2 x 2 Lax matrix L(A) = K + T(X)K_T~ 1 (-X) gj. The tensor P x (1631) 
differs from the Poisson tensor for the Kowalevski gyrostat, which appears from the linear 
r-matrix algebra 



The 2x2 Lax matrix for the Goryachev-Chaplygin gyrostat looks like [8] 





T(A) 



2 





_2 

e 2 



X + 2J 3 +p a 
a 



T(A) 



e 2 






2 



(6.4) 



Here p, q are additional dynamical variables, a is an arbitrary number and T(X) is given by 

(EU). 

Substituting this matrix into the Sklyanin bracket (11. ip and into the brackets (I4.7p - (I4.8I) 
one gets the compatible Poisson tensors on the extended phase space e*(3) x (p,q) 



( 



Pn 



W 
Go 



Po 



\ 



2i 
* J 



(6.5) 



and 



/ 



Pi 



Pi W x 



Pi 



\ 



-2x 3 J 2 + 2px 2 + 8x 2 J 3 
2x 3 Ji — 2px\ — 8xi J 3 
2x\J 2 — 2x 2 J\ 
2J 2 (p + 3J 3 ) 
—2ax 3 — 2pJ\ — 6Ji J 2 
2ax 2 



-2x 2 
2x x 


-2 J 2 + ix 3 e q 
2 Ji - x 3 e 9 
-i(xi + ix 2 )e q 



\ 



2i((xi - ix 2 )e 9 - 2 J 3 - p - ae«) 
* 



/ 



(6.6) 



which satisfy equations (12. ip at C 2 = only. Here P and P x are given by (16.21) and (16. 3p . 
The Hamiltonians Hi from the trP(A) = A 3 + PiA + X 2 H 2 + H 3 are 



^i 
P 2 

#3 



-(J x + J 2 + 4J 3 2 + 2pJ 3 - 2ax x ) 
-{2J 3 +p){J 2 l +Jl)-2ax 3 J l . 



The obtained tensors Po and P x are degenerate and, therefore, the Hamiltonians Hi reproduce 
the Frobenius chain in the following form 



P 1 dH l = P {dH i+x + Ci dH x ) 



1,2,3, 



(6.7) 
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where H A = and q are coefficients of the polynomial Ajv(A) = A(X) + B(X) — C(A) — D(X) 

At p = p and q = matrices Go (16.51) and G\ (16. 6p are (generically) non-degenerate. So, 
the Dirac procedure can reduce pencil Pq + A Pi to a new Poisson pencil P D + APj° on e*(3) 
defined by 

P k ° = Pk+ (W k G~ k l Wl) p=p q=0 , k = 0, 1. 

Here Pq = Pq is canonical Poisson tensor ( 16.21) and Pi is given by ( 16. 3j) . This reduction 
procedure preserves equations ( 16.71) for the reduced integrals of motion. 

7 Conclusion 

We present a family of compatible Poisson brackets (I2.10p . that includes the Sklyanin bracket, 
and prove that the Sklyanin variables are dual to the special Darboux-Nijenhuis coordinates 
associated with these brackets. The application of the r-matrix formalism is extremely useful 
here resulting in drastic reduction of the calculations for a whole set of integrable systems. 

The construction can be generalized to other r-matrix algebras. Remind, if one substi- 
tutes T(A) = 1 + sL(X) + 0(e 2 ), r = er into ([11]) and let e -> one gets a linear bracket. 
Then if T(A) satisfy the Sklyanin bracket (11.11) . then the matrix T(A) = T{X)K^T^ 1 {— A) 
obeys to the reflection equation algebra [10] • The corresponding compatible brackets for the 
open generalized Toda lattices was considered in [15] . 

Moreover, the whole construction can immediately be transferred to the quantum case 
because r-matrices in (12.101) became dynamical matrices at k > 1 only. 

We would like to thank I.V. Komarov and V.I. Inozemtsev for useful and interesting 
discussions. The research was partially supported by the RFBR grant 06-01-00140. 
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